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We analyze random unitary evolution of a qubit within memory kernel approach. We provide 
sufficient conditions which guarantee that the corresponding memory kernel generates physically 
legitimate quantum evolution. Interestingly, we are able to recover several well-known examples 
and to generate new classes of nontrivial qubit evolution. Surprisingly, it turns out that a class of 
quantum evolutions with memory kernel generated by our approach gives rise to the vanishing of a 
non-Markovianity measure based on the distinguishability of quantum states. 
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I. INTRODUCTION 

Dynamics of open quantum systems plays an impor¬ 
tant role in the analysis of various phenomena like dis¬ 
sipation, decohcrence and dephasing 01- The usual 
approach to the dynamics of an open quantum system 
consists of applying the Born-Markov approximation 0 
which leads to a local master equation for the Markovian 
semigroup 

Pt = L[p t ] , (1) 

where pt is the density matrix of the investigated system 
and L is the time-independent generator of the dynamical 
semigroup defined as follows 

L[p] = -i[H, p] + ^ X ([^> P V «] + [ V "P> V «}) ■ ( 2 ) 

a 

Here H denotes the effective system Hamiltonian, and 
V a represent noise operators 0 . 0 . We call Q the 
GKSL form (Gorini-Kossakowski-Sudarshan-Lindblad). 
The solution of (JT]) defines the Markovian semigroup 

p t = A t [p] = e tL p , (3) 

where p is an initial state. The dynamical map A t = e tL 
is completely positive and trace-preserving (CPTP) 0,Il¬ 
ls] . The Born-Markov approximation assumes weak inter¬ 
action and a separation of time scales between the sys¬ 
tem and its environment. Such approach works perfectly 
well for many quantum optical systems 0]sj. When the 
above assumption is no longer valid the description based 
on ED is not satisfactory. Recent technological progress 
and modern laboratory techniques call for a more refined 
approach which takes into account memory effects com¬ 
pletely neglected in the description based on Markovian 
semigroups. In recent years we observed an intense re¬ 
search activity in the field of non-Markovian quantum 
evolution (see the recent review 0, a collection of arti¬ 
cles in 0 and a recent comparative analysis in 0)- 

There are basically two approaches which generalize 
the standard Markovian master equation (ED: time-local 
approach replaces L by a time-dependent generator L t . 
Interestingly, if for all t the time-dependent generator has 


the standard GKSL form ©, then A t = T exp( L u du ) 
defines the so-called divisible dynamical map jl2l . [131 ] 
which is often considered as the generalization of Marko- 
vianity (see 01] for a generalization of the notion of di¬ 
visibility). The second approach is based on the nonlocal 
Nakajima-Zwanzig equation 0 (see also 0 ) 

Pt= Kt-rPrdT, ( 4 ) 

Jo 

in which quantum memory effects are taken into account 
through the introduction of a memory kernel K t . This 
means that the rate of change of the state pt at time t 
depends on its history (starting at t = 0 ). The Marko¬ 
vian master Eq. (1) is reobtained when K t = 2S(t)L. 
The time-dependent kernel is usually referred to as the 
generator of the non-Markovian master equation. Equa¬ 
tion ED applies to a variety of situations (see eg. 07|). 
Because of the convolution structure of ED the time-local 
approach is often called time-convolutionless 0130 . 
The structure and the properties of ED were carefully 
analyzed in [20l429j | . In particular the generalization of 
Markovian evolution to the so-called semi-Markov was 
investigated within the memory kernel approach by Bu- 
di ni Ell and Breuer and Vacchini [23} (see also discussion 
in 

In a present article we study random unitary evolution 
of a qubit within the memory kernel approach. In partic¬ 
ular we address the following problem: what is the struc¬ 
ture of the corresponding memory kernel Kt which leads 
to the legitimate CPTP dynamical map A t . The article 
has the following structure: in Section II we recall basic 
facts about random unitary evolutions and in Section III 
we formulate the sufficient condition for K t to guarantee 
legitimate physical evolutions. In Section HVl we examine 
the issue of Markovianity. Surprisingly, it turns out that 
a subclass of quantum evolutions with memory kernel 
generated by our approach gives rise to the vanishing of 
a non-Markovianity measure based on the distinguisha¬ 
bility of quantum states [3] , Section V illustrates our 
approach with several examples. Final conclusions are 
collected in Section VI. 


2 


II. RANDOM UNITARY QUBIT EVOLUTION 

A quantum channel £ : 13(H) -A B(H) is called random 
unitary j32j if its Kraus representation is given by 

£[X] = Y,p^ XU l > ( 5 ) 

k 

where Uk is a collection of unitary operators and {pk} 
stands for a probability distribution. In this article we 
consider a random unitary dynamical map A t defined by 


3 


A t [p\ = > ^Pa(t)cr a pa- a , 


ct—0 


( 6 ) 


where a a are Pauli matrices with <tq = I 2 [331 ]. Initial 
condition A t= o = 1 implies p Q (0) = S a o- Recently a time- 
local description based on the following master equation 
was analyzed Him 

A* = L t A t , (7) 

where Lt is a time-local generator defined by 
3 

Lt[p] = 7 k(t) ((TkpcTk - p), ( 8 ) 

fc=1 

with time-dependent decoherence rates 7 k(t). One asks 
the following question: what are the conditions for 7 k(t) 
which guarantee that the solution At = exp (J Q L T dr) 
provides a legitimate dynamical map? Note, that the 
solution defines a random unitary evolution with p a (t) 
given by 

3 


Pa(t) — , y ] H a pXp(t) 


(9) 


3=0 


where H a p is the Hadamard matrix 


H = 



( 10 ) 


III. CONSTRUCTION OF LEGITIMATE 
MEMORY KERNELS 

In this article we analyze the nonlocal description 
based on the following memory kernel equation 


with 


At — J K t - T h. T &T , 

0 


K t[p\ = kl ^> ( a iP ai _ P) > 


(15) 


(16) 


where ki(t) (i = 1,2,3) represent nontrivial memory ef¬ 
fects. Note, that equation m considerably simplifies 
after performing the Laplace transform 


A., = 


1 


s -K s 


(17) 


where A s := e~ st A t dt and similarly for K s . The ques¬ 
tion we address is: what are the conditions for ki(t) which 
guarantee that the solution A t provides a legitimate dy¬ 
namical map1 

Denoting by n a (t) the eigenvalues of K t: 

K t [a a ] = K a (t)cr a , (18) 

equation (fl5l) gives rise to the following set of equations: 

A i(t)= [ Ki(t - r)A»(r)dr, * = 1,2,3. (19) 

Jo 

Note, that Ko(t) = 0 and hence Ao(t) = 1 = const. In 
terms of the Laplace transforms At(s) and k;(s) one finds 

1 


A,(s) = 


s - ki(s) 


( 20 ) 


In terms of Ai(s) conditions (HI (O may be equiva¬ 
lently reformulated as follows: 


and Xp(t) are time-dependent eigenvalues of A t 
Ai[(Ja] — X a (t)(T a , 

(11) 

—b Ai(s) + Ai(s) + A 2 (s) 
s 

is CM 

which read as follows: Ao(i) = 1 and 

and 



Ai (t) = exp(—2[T 2 (f) + T 3 (t)]), 


—b X 3 (s) — Al(s) — A 2 (s) 

is CM 

X 2 (t) = exp(—2[Ti(t) + r 3 (t)]), 

(12) 

s 

1 - 


X 3 (t) = exp(-2[Ti(t) + T 2 (t)]), 


- + A 2 (s) - Ai(s) - A 3 (s) 
s 

is CM 

with Tfc(t) = J*jk(T)dT. Now, the map © is CP iff 
p a (t) > 0, which is equivalent to the following set of 

—b Ai(s) — As(s) — A 2 (s) 
s 

is CM 


conditions for As [34|, [35 


1 + Ai (t) + X 2 (t) + A 3 (f) >0 , 


where CM stands for a completely monotone function 
(13) [37]], i.e. a smooth function / : [0,oo) -A R satisfying 

the condition 


Ai(t) + X 2 (t) < 1 -I- A 3 (t), 
X 3 (t) + Xi(t) < 1 + A 2 (t), 
X 2 (t) + X 3 (t) < 1 + Xi(t). 


( 14 ) 


(_ 1 ) "d^ /(s) -°’ S -°’ n = 0 ’ 1 ’ 2 ’--- ( 23 ) 

The equivalence of (THl) and results from the follow¬ 
ing 


and 
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Theorem 1 (Bernstein’s Theorem) A function 
/: [ 0 , oo) —>• K is completely monotone on [ 0 , oo) if and 
only if it is a Laplace transform of a finite non-negative 
Borel measure p on [0, oo), i.e. f is of the form 

POO 

f(s) = / e~ st d p(t). (24) 

Jo 

Note that the initial condition po(0) = 1 and Pk(f>) = 0 
for /c = 1,2,3 is equivalent to Afc(0) = 1 due to 

Theorem 2 (Initial Value Theorem) Let f(s) be the 
Laplace transform of f(t). Then the following relation is 
true: 


1 

—b 

1 

> 

1 

a 1 

02 


03 

1 

-b 

1 

> 

1 

02 

03 


Oi 

1 

— + 

1 

> 

1 

03 

ai 


02 


define a legitimate memory kernel 
1 3 

K s [p\ = -^K a (s)<Tfctr[a- fc p] , (31) 

k =1 


lim /(f) = lim sf(s) (25) 

t—tO s—yoo 

it is equivalent to 

lim sAfe(s) = 1 , (26) 

s—»■ OO 

for k = 1,2,3. This way we have proved 

Theorem 3 The map A s represented by the following 
spectral decomposition 

1 3 

A s [p] = - ^2 \ a {s)cj a tr[cj a p\ , (27) 

a=0 


with Xq (s) = 1/s, defines the Laplace transform of a le¬ 
gitimate map A t if and only if conditions \21\) . H22\) and 
11?61) are satisfied. 


It is worth emphasising that there are few analytical tools 
for dealing with CM functions, which is due to the fact 
that an infinite set of conditions (l23l) must be verified. 
Nevertheless, we found an important class of CM func¬ 
tions giving rise to CPTP dynamics with a straightfor¬ 
ward interpretation. To present them, let us first ob¬ 
serve that CM functions have the following two proper¬ 
ties, which will not be proved: 

Property 1 Let f and g be arbitrary completely mono¬ 
tone functions. Then 

1. f ■ g is CM, 


2. af + fig is CM for any a, /3 > 0, 

Property 2 If so > 0 then is CM. 

We are now ready to prove our main result: 

Theorem 4 Let W (s) be a function such that ^ 
CM. Then the functions 


Kk(s) = — 


a k W(s ) - 1’ 
with ai, 02,03 > 0 such that 


k = 1,2,3, 


1 

s 



1 1 

W(s) |_Ql 0,2 



is CM , 


(28) 


(29) 


i.e. the corresponding X &(s) satisfy \21\ ) 
Proof: note that formula (l28l) implies 

1 


and 


X fc (s) = - ( 1 - 


and hence 


a k W{s) 


—b A 3 (s) — Ai(s) — A 2 (s) 
s 

11/1 1 1 

s W(s) \ai a 2 03 , 


(32) 


(33) 


which proves that | + As(s) — Ai(s) — A 2 (s) is CM due 
to the fact that is CM. Similarly one proves the 

remaining conditions (11411 . □ 

Note, that since - . is CM, hence, due to the Bern- 
stein theorem, it is the Laplace transform of a positive 
function. Hence 


W(s) = 


1 


/OO ’ 


(34) 


where /(s) is the Laplace transform of /(f) satisfying 
fo f( T )dr > 0 for all t > 0. One finds 


Kfc(s) = 

Note, that condition 


1 

ai 


1 

a 2 


~ s f(s) 

Ok - f(s) 
implies 
1 

03 


j J f( T )dr < 4 . 


(35) 


(36) 


Hence to summarize: our class is characterized by a single 
function /(f) and three numbers ai,a 2 ,a 3 > 0 such that 
F(t) = fg f{r)dr > 0 and conditions (13(1 and (1551) hold. 
One finds for p a (t): 


Pi(t) = 


1 


P2(t) = Z 
P3(t) = l 


1 

02 

1 

03 

1 

ai 


1 

03 

1 

0\ 

1 

02 


1 

ai 

1 

02 

1 

03 


F(t) 
F(t ) 

m 


( 37 ) 
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and po(t) = 1 —p\(t) —p 2 {t) — p 3 (t). In particular, taking where Xk '■= tr [p<Jk\ and the relaxation times are defined 
ai = a ,2 = a and a 3 = oo one finds via 


ki(s) = k 2 (s) = —, n 3 (s) = 0 , (38) 

a~ f{s) 

and hence 

Ms) = k(s) = 0 , k 3 (s) = , (39) 

2 a- f{s) 

gives rise to the legitimate memory kernel 

K t \p\ = k 3 (t)(a 3 pa 3 - p), (40) 

with arbitrary f(t) and a > 0 satisfying additional con¬ 
dition 


0 < F(t) := f f{r)dT < 2a, (41) 

Jo 

for allt > 0. The corresponding solution reads 

Po{t) = 1 - Tir(t) > 

Piit) = P2(t) = 0 , (42) 

PM = ±F<f) . 

This approach resembles very much the semi-Markov 
construction [ 2 ^, for any f(t) > 0 satisfying 

f 0 f(t)dt < 1 the memory kernel (T40l) with 


Ti = , T 2 = , T 3 = . (48) 

72 +73 73 + 71 71 + 72 

It is well known Q that complete positivity is equivalent 
to the following set of conditions upon T^: 


1 

+ 

1 

> 

1 

T~i 


% 


n 

1 

+ 

1 

> 

1 

% 


T3 


t[ 

1 

-f 

1 

> 

1 

T3 




T2 


(49) 


It is therefore clear that condition (1301) is an analogue 
of ([49]) . Note that condition (l30l) means that there exist 
bi,b 2 ,b 3 > 0 such that 


1 1 _ 1 1 

2 a 3 62 b 3 

1 1 _ 1 1 

2 a 2 63 + b \ 

1 1 _ 1 1 

2 a 3 bi 62 


(50) 


Now, it terms of 61 , 62 , b 3 our result may be reformulated 
as follows 


Corollary 1 For any b±,b 2 ,b 3 > 0 and the function f(t) 
satisfying 


k 3 {s ) 


s f(s) 

i-/(T 


(43) 


0 < F(t) := 


[ f{r)dr < 
Jo 



1 

^2 



(51) 


gives rise to CPTP evolution. In this case one finds 


and 


Po(t) = ^[1 + Ai(t)] , 

Pi(t) = p 2 {t) = 0 , (44) 

P3(t) = \\ l ~ Ai(*)] , 

where 

A 1 ( S ) = A 2 ( S ) = ^ 4 ± T- ( 45 ) 

/(s) - 1 


It is therefore clear that our approach goes beyond the 
semi-Markov construction. 

Let us recall that Markovian semigroup generated by 


L[p] 


1 

2 


3 

'Y^ikiakpvk - p\, 

k=l 


the corresponding Bloch equation reads 


Xkit) 



(46) 


(47) 


lim /(s) = 0 , 

s—^oo 

the memory kernel defined by 

sf(s) 


Kfe(s) = — 


(52) 


(53) 


afc - f(s) 

defines legitimate quantum evolution. Moreover one has 


Pkit) = ji-F(t) , 


(54) 


andp 0 (l) = l-pi{t) -p 2 (t) -p 3 (t). 


Let us observe that it is very hard, in general, to invert 
formula (f35l) to the time domain. Now, we provide a 
family of W (s) which enables one to easily compute «q(f) 
and have the memory kernel in time domain. 

Theorem 5 Let W{s) be a polynomial 

kP(s) = (s + zi)... (s + z n ), (55) 










5 


with Zi > 0. //ai, < 22,02 satisfy \3(A) and 

(l 1 1 

A 1-1-1- 

4 \ ai a 2 a 3 


” 1 

n^z 


i=1 


f/ien Ki(t) defined via 
kernel. 


(56) 

define a legitimate memory 


Proof: It is clear that it is enough to prove m- 
Lemma 1 One has the following decomposition 

n i —1 \ 

7 = 1 \ 


1 

sYVUs + z*) 

where 


= A 


--Y 

s 1= i ni= 1 (a + *i) > 


A = 


nr= 


(57) 


(58) 


i=l ** 


For the proof see the Appendix. Now we show that con¬ 
dition (Oil) holds. According to (1571) one has 


+ Ai(s) + A 2 (s) + A 3 (s) 


1 A 

1 1 1 

= - 4- 

— + - 1- 

s V 

ai a 2 a 3 _ 

1 (a 

' 1 1 1 ' 

= - 4- 

- \ - 1 - 

s V 

ai a 2 a 3 _ 


1 


W{s) 

1 


(59) 


1 

ai 


1 

02 


1 

o 3 


E 


nti - 


nili ZifriYlhis + ZiY 


A second term in (1591) is CM due to the fact that it is a 
sum of CM functions. Hence, if condition (l56l) is satisfied 
then (1711) holds. □ 

Note, that 


tii{s) = — 


1 


CLkW (s) - 1 Ofc (s - Si) . . . (s - Sm) ’ 

(60) 

where {si,...,s m } are the roots of the polynomial 
(afcW(s) — 1). It is therefore clear that formula (1551) may 
be easily inverted to the time domain. 

Remark 1 Note, that W (s) defined in 1,551) implies that 


W(s 


is CM and hence - 


W(s 


is CM as well. 


IV. CHECKING FOR NON-MARKOVIANITY 


Let us recall that according to [3l| the evolution rep¬ 
resented by A t is non-Markovian if the condition 


dt 


||A t [pi - P 2 ] | |tr < 0 , 


( 61 ) 


is violated for some initial states pi and p 2 - One defines 
31| a well-known non-Markovianity measure 


/ d 

— ||A t [pi - p 2 ] 11tr dt , 
dt 


(62) 


where the integral is evaluated over the region where 
^||A t [p! — p 2 ]||tr > 0. Now, it has been proved H3| 
that for random unitary qubit evolution if all eigenvalues 
A kit) > 0, then (1CT1) is equivalent to 

^Afc(f) < 0 ; = 1,2,3. (63) 

Proposition 1 For ai,a 2 ,a 3 satisfying 130 1 ) and 
W(s) = Yy where f(s) is CM and 


/(r)dr < 

^min? ^min — min{ai,a 2 ,a 3 }, (64) 

the corresponding memory kernel gives rise to the dynam¬ 
ical map A t such that A/BLp[At] = 0. 

Proof: Let us observe that condition (l64l) implies (l29l) . 
Indeed, from (1671) one has 



— + — + — ) [ f(r)dr < 3, (65) 

ai a 2 a 3 y J 0 


and hence the condition 


Afc(t) = 1 - — [ /(r)dr > 0 
a k Jo 


follows. Now, observe that 

t 


due to (1571) . Hence it is sufficient to show that ^A k(t) < 
0. It is clear ^A k(t) < 0 if and only if 1 — s\k{s) is CM 
and hence taking into account (12U1) it is equivalent to the 
requirement that — kfc(s)Afc(s) is CM. One has therefore 


- K k (s)Xk{s) = ( 66 ) 

^k 

which ends the proof since f(s) is CM and ak > 0 □ 

Remark 2 It W(s) = (s + z\)... (s + z n ) with Zk > 0 
and ai,a 2 ,a 3 satisfying 13(A) together with 


n 

T*>-, fc = 1,2,3, (67) 

i=i ak 

then the corresponding dynamical map A t satisfies 
A/BLp[A t ] = 0. 

Remark 3 It was shown f!2 . \sA] that BLP condition 
m is equivalent to so-called P-divisibility. This means 
that 


A t = V t , s A s , (68) 

and for any t > s the propagator Vt :S is positive (but not 
necessarily completely positive). 

Interestingly, our construction provides a class of legiti¬ 
mate random unitary qubit evolutions generated by the 
nontrivial memory kernel but still satisfying BLP condi¬ 
tion (1511) . (cf. also [30j). It is clear that to violate dSU 
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one needs a more refined construction such that is 
not CM but is already CM. It deserves further 

analysis. 

Consider now the question of CP-divisibility which is 
fully controlled by the local decoherence rates in ([ 8 ]). One 
may easily compute them in terms of f(t ): 


mo-f' - 1 


1 


1 


^ M fit) 

72 (i) = —r 


4 \ai — F(t) a 2 -F(t) a 3 -F(t)J' 

1 1 1 


73 (f) 


1 \ai — F(t) a 2 -F{t) a 3 -F(t)J ’ 
f{t) ( 1 1 1 


4 \ai — F(t) a 2 — F(t) a 3 — F(t)J 


The dynamical map A* is CP-divisible iff 7 k(t) > 0 for 
k = 1.2,3. Let us assume that 


a l ^ a 2 a 3- 


(69) 


Proposition 2 If a\,a 2 ,a 3 and /(f) > 0 satisfy condi¬ 
tions $29 J) and i t.'?dll f/ie corresponding memory kernel 


Kk(s) = -- 


s /0) 


(70) 


a fe - /(s) 

leads to a CP-divisible dynamical map iff 

F(t) < ai - \/(a 2 - Oi)(a 3 - ai). (71) 

Proof: Due to (15^1) it is sufficient to show that 71 (f) > 0 
which, for /(f) > 0 , is equivalent to 


-1 


1 


1 


ai - F(t) a 2 - F(t) a 3 - F(t) 


> 0. (72) 


Let us assume that F(t) < ai, which means, that 71 (f) 
is not singular. Inequality (1771) is satisfied iff 


with 


F(t) € (—oo,F_] U [F+,+oo) 


F± = ai ± V (a 2 - ai)(a 3 - ai). 


Now, taking into account that F(t) < a\ one finally 
proves ED- □ 

This Proposition shows that positivity of the function 
/(f) is not sufficient for CP-divisibility. One needs an 
extra condition CD which involves not only /(f) but 
{ai, 02, <23} as well. 


V. EXAMPLES 

Example 1 Consider the simplest case with a polyno¬ 
mial of degree one 


(73) 


with z > 0. One finds 

fik(s) = - 


a k (s + z) - 1 ’ 
and the inverse Laplace transform gives 


Kfc(f) = -- ( 5(t) - 


z — 




-\z~-±]t 
3 L a k 1 


(74) 


(75) 


Note, that if ak = 1 /z, then the dynamics is purely local. 
One easily finds 


1 

za k 


A fc (t) = 1-(1 — e~ zt ), 


(76) 


and finally the solution for pkit) is defined by © with 


F{t) = 7(1 — e~ zt ) . 

z 


(77) 


Note, that condition (EDI implies the following relation 
between z and ai,a 2l a 3 : 


. 1 1 1 

4 z >-b-t-, 

®1 02 ^3 


(78) 


which guarantees that po(t) >0. In the symmetric case 
ai — a 2 = a 3 — a one finds pi(t) = p 2 (t) = p 3 (t) p(t) 
with 




(79) 


and Po(t) = 1 — 3p(f) with 4za > 3. One finds that 
asymptotically 


"» (i) -* 1 - TTa 


(80) 


Note that for za > 1 one has asymptotically po(oo) < 
1/4. This property cannot be reproduced within the local 
approach with regular generators L t . Indeed, it follows 
from m (see also W2 l for more details) that 


Po{t) — j[l + Ai (f) + A 2 (f) + A 3 (f)], (81) 


and hence, using one finds 


Po{t) > j- 


(82) 


This example shows that local and memory kernel ap¬ 
proaches may lead to essentially different evolutions. 

Example 2 Consider now the same polynomial W ( s ) = 
s + z but let z = 2c > 0. Moreover 


1 


ai — a 2 — 


a 3 = 


1 

2c 


(83) 


One finds 


ki(s) = h 2 {s) = -■ 


sc 


W(s) = s + z, 


s + c 


k 3 (s) = - 2 c, 

























7 


and hence 


Example 4 Let 


Ki(t) = « 2 (t) = —cS(t) + c 2 e ct , k 3 (i) = —2 cS(t). 


W(s) = s 2 +u 2 . 


(92) 


Finally, one finds the following formula for the memory 
kernel 

Q 

K t[p\ = -<5(t)[aipoi + a 2 pa 2 - 2p] 

- ye _rt [cr 3 pcr 3 -p], (84) 

One ftas 

Ai(t) = A 2 (i) = i (1 + e“ 2ct ) , Aa(t) = e- 2ct . 


Note i/iai i WXsJ since 

11 _ 1 1 ( w A 

sTE(s) ws\s 2 +u; 2 / 

is the Laplace transform of f* sin(wr)dr which is positive 
for all t > 0. Condition 1291) implies 

1 1 1 9 

— + — + — < 2w 2 . 93 

ai a 2 a 3 


Interestingly, this evolution reproduces time-local descrip¬ 
tion with 

7i (*) = 72 (i) = 73 (i) = tanh(ct). (85) 

as discussed in fOdi l. It was shown f§di] that A t is a convex 
combination of two Markovian semigroups and A[ 2 ^ 
generated by 

Q 

L k[p\ = -^WkP^k - p\ ; k= 1,2, (86) 

</iai is, 

A* = i (e tLl + e* i2 ) . (87) 

This simple example shows that a convex combination 
of Markovian semigroups leads to a quantum evolution 
displaying essential memory effects. 


Example 3 Consider now a polynomial of degree two 


W(s) = (s + ci)(s + c 2 ), (88) 


with c 2 > ci > 0. Our construction gives rise to a legiti¬ 
mate memory kernel if condition \30\) holds and 


1 

4 cic 2 >- 1 - 

a 1 


1 

a 2 


1 

a 3 ‘ 


One finds 


Kfc(s) 


1 s 

afc (s + ci)(s + c 2 ) - i 
1 s 

ak (s + si)(s + s 2 ) ’ 


(89) 


(90) 


with 


1 

Sl + s 2 = Cl + C 2 , SiS 2 = Cic 2 -. 

ak 

Hence the solution has the form with the function 
F(t) given by 


F(t) 


1 

c 2 - Ci 




— e 


-C- 2 .t 


(91) 


The corresponding eigenvalues of the memory kernel read 
nfit) = - — cos ( Ju) 2 - — t J , (94) 

for w 2 > 1 /ai, and 

Ki{t) = — — cosh ( J —— w 2 t ) , (95) 

at VV J 

for uj 2 < 1 /ai. Moreover one finds 

A k (t) = H-^[cos(wi) - 1] , (96) 

a fe w 2 

and hence 

P i 0 ) = ( ~ + --) [1 - cos(wi)] , 

4ar \a 2 a 3 ai J 

P 2 {t) = (- + - --) [1 - cos(wt)] , (97) 

4ar \ a 3 ai a 2 J 

Pa(t) = T~2 {— + ---) I 1 - cos M)] , 

4ar \a 1 a 2 03 / 

together with po(t) = 1— pi(t) — p 2 {t)— p 3 (t). In particular 
taking 

1 

ai=a 2 = —, a 3 = 00 , (98) 

one finds 

Ki(t) = K 2 (t) = -w 2 , « 3 (t) = 0 , (99) 

and hence 

ki{t) = k 2 {t) = 0 , n 3 (t) =, ( 100 ) 

which proves that the constant (time independent) 

K t[p\ = ^(cr 3 pa 3 ~ p), ( 101 ) 

provides a legitimate memory kernel for arbitrary k = 
uj 2 > 0. Moreover one finds for the local decoherence 
rates 











= ui sin(o;f) /_ -1 _ 

4 Vaiw 2 — 1 + cos(wt) 

_ 1 _ + _1_\ 

a 2 W 2 — 1 + cos(wf) a^uj 2 — 1 + cos (cot) J ’ 

and similarly for 72 (t) and 73 (f). Vote, that if for some k 
one has akbJ 2 < 1 £/ien local decoherence rates are singu¬ 
lar and hence in this case the non-local approach is more 
suitable. 

VI. CONCLUSIONS 

We analyzed random unitary evolution of a qubit 
within memory kernel approach. Our main result formu¬ 
lated in Theorem [5] allows to construct legitimate mem¬ 
ory kernels leading to CPTP dynamical maps. The power 
of this method is based on the fact that 1 ) it allows to 
reconstruct well known examples of legitimate qubit evo¬ 
lution, 2 ) the structure of polynomials R 4 (s) enables one 
to perform the inverse Laplace transform and to find a 
formula for the kernel in the time domain. The math¬ 
ematical analysis heavily uses the notion of completely 
monotone functions. These functions are not commonly 
used in theoretical physics and knowledge of their prop¬ 
erties is rather limited. There are no known effective 
methods allowing to check whether a given function is 


CM. We stress that Theorem [4] provides only a sufficient 
condition and further analysis is needed to cover physi¬ 
cally interesting cases which do not fit the assumptions of 
the Theorem. Interestingly, it turns out that the quan¬ 
tum evolution with a memory kernel generated by our 
approach gives rise to vanishing non-Markovianity mea¬ 
sure based on the distinguishability of quantum states 
[3lj . We also have shown when the corresponding dy¬ 
namical map is CP-divisible. It shows that the evolution 
satisfying nonlocal master equation does not necessarily 
lead to a non-Markovian evolution. It would be also in¬ 
teresting to analyze the relation between semi-Markov 
evolution and the one governed by our approach in more 
detail. 
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Appendix: proof of Lemma [T] 

Let us observe that (1571) may be represented in the 
following form 


• s n"=i(s + ^) 


= A 


n"=i( s +Zi) ~ s (n;u(s ++~i nr=3( s ++ • • ■ nj=i ~?( s + z n )+ riy=i 

«nr=i (*+*) 


( 102 ) 


therefore, to prove the Lemma it suffices to show that 


J 


n* 


n 

n(*+*o-s 

i— 1 


+ Zi) + Z\ + Zi) 


2 


i —3 


n —2 n —1 \ 

• n z i ( s + z «)+n z i 

i=i 1=1 / 


We will prove this by induction. For n = 1 it is clear that while RHS reads 
LHS=RHS=,ti. We assume that (11031) is true for n and 
prove it is also true for (n + 1). LHS may be written as 

n 

LHS - 11 Zi • z n ~ |-i, 

i= 1 


(103) 


(104) 
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RHS = 


[[(* + z i X s + z n+ 1) — J^(s + z i)(s + Z n + 1 ) + Z\ J^[(s + z i)(s + Z n + 1 ) + . . . 

i— 1 i —2 i=3 

n— 2 n— 1 n— 1 

+ z j(s + Z n )(s + z n - i-i) + Zj(s + z n + 1) + Zj ■ z n ^ = 

1=1 1=1 1=1 

/ n / n n n—2 

= (s + Z„_|_i) I J^[(s + ^i) — s I JJ(s + Zi) + Z 1 JJ(s + Zj) + . . . Zj(s + Z„) + J 


n —1 


n— 1 


L 2=1 


n 


2=2 


2=3 


3 =1 


J=1 


n 


n 


- s JJ ZjZ n = s Zi + z„+i ]^[ Zi - s ]^[ nZi, 
1=1 i=1 


(105) 


i=l i=l 


r 


which proves that RHS=LHS. 
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